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The Langevin equation with a multiplicative Levy white noise is solved. The noise amplitude 
and the drift coefficient have a power-law form. A validity of ordinary rules of the calculus for the 
Stratonovich interpretation is discussed. The solution has the algebraic asymptotic form and the 
variance may assume a finite value for the case of the Stratonovich interpretation. The problem of 
escaping from a potential well is analysed numerically; predictions of different interpretations of the 
£SJ ' stochastic integral are compared. 
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I. INTRODUCTION 



The Langevin equation can not always be expressed by means of a deterministic drift term, supplemented by a time- 
dependent stochastic force (an additive noise). A physical quantity, which is represented by the random component, 
may require the noise to depend on the stochastic variable itself. In the Langevin description, that dependence 
emerges as a variable noise amplitude (a multiplicative noise). The multiplicative noise emerges also in descriptions 
of complicated systems, as a result of the elimination of fast degrees of freedom. The stochastic equation is then of 
I ' the form 



x = F(x) + G(x)ri(t), (1) 



a ■ 
•i— > . 

where F(x) and G(x) are given functions. The stochastic force r/(t) is uncorrelated, {r](t) r l(t')) = S(t — t'), and it is 
characterised by a given probability distribution. In the present paper we assume that 77 has the symmetric stable 
^p 1 ' Levy distribution defined by the Fourier transform 

-a' 

L a (k)=cxp(-\k/a\ a ), (2) 

O 

where a is the order parameter and a scales the distribution. In Sees. II and III we assume a = 1. The case a = 2 
corresponds to the normal distribution which is well known in the context of multiplicative processes [l|, ■ 
£N) ■ The general and stable Levy processes, for a ^ 2, exhibit long tails of the distribution which makes the moments 
divergent. They are frequently encountered in nature, since long jumps are associated with a complex structure of the 
environment, in particular with long-range correlations. Examples can be found in biological physics [3|], disordered 
media [|| and finance 043 • A master equation description of thermal activation of particles within the folded polymers 
[H also involves the multiplicative Levy noise in a sense that the equation is fractional (Levy jumps) and it contains 
a variable diffusion coefficient. However, a direct Langevin representation of the topological complexity problem is 
unknown Q. Since the complex environment is usually nonhomogeneous, one can expect that the Levy noise in 
the Langevin equation is rather multiplicative than additive. Therefore formalisms, which are supposed to describe 
complex processes and which do that in terms of the additive noise alone, may miss essential features of the problem. 
For example, in the field of finance, the standard Black-Scholes equation contains the additive Gaussian noise. Eq.([T]) 
in its general form could be an important generalisation of that equation 0, |(| . 

Eq.(Q} is not sufficiently defined for the white noise because it is not clear at which time G(x(t)) should be evaluated. 
In the following, we define the stochastic integrals, connected with Eq.([T]), as Riemann integrals. According to 
Stratonovich, one assumes 



/ G[x(r)}d V (r) =£G 



x(U-i) + x(U) 



[V(U) - v(U-i)l (3) 



where t = rn and r = tj — is a time step. This interpretation is appropriate for many physical phenomena since 
it constitutes a white noise limit for correlated processes. In this case Eq.ffl c an be solved like usual differential 
equation, which can be rigorously proved if 77 has the convergent variance [10. Illf. In particular, one can introduce a 
transformation 
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which leads to the Langevin equation with the additive noise: 

y = F(y)+ V (t). (5) 

Alternatively, we can simply assume 

/ G[x(r)}dT,(r)=J2G[<ti-i)MU)-v(U-i)], (6) 
J ° i=i 

which formula defines the Ito interpretation. Predictions of Eq.© in both interpretations are different but in the case 
a = 2 in Eq. © there is a simple relation between them: the difference resolves itself to the spurious drift [l(| • For a < 2 
such a relation does not exist. The Stratonovich interpretation predicts a dependence of the probability distribution 
on the noise amplitude which may change the diffusion properties of the system; in particular the accelerated diffusion, 
in the case of the additive noise, can change to the sub-diffusion. That problem is discussed in Ref.[l2[ for the case 
without drift and for the linear drift. In the Ito interpretation, in turn, shape of the distribution tail is not affected 
by the amplitude [la ]. 

In this paper we discuss properties of the Langevin equation which is driven by the multiplicative Levy noise and 
nonlinear forces, in particular the problem of escaping from a potential well. In Sec. II, properties of stochastic integrals 
for the stable Levy processes and those with truncated distributions are compared. The Fokker-Planck equation for 
the problem of an algebraic, nonlinear potential is solved in Sec. III. The escape from the potential well, understood 
as the first passage time problem, is calculated in Sec. IV and results for both interpretations of the stochastic integral 
are compared. 



II. STABLE LEVY DISTRIBUTIONS VERSUS TRUNCATED ONES 

A well known property of the Stratonovich integral © allows us to apply standard rules of the calculus and then to 
reduce Eq.© to an equivalent equation with the additive noise. It can be proved [Hi [HI for the normally distributed 
noise, i.e. on the assumption that increments are independent and the variance is finite. For the general Levy stable 
processes the latter condition is not satisfied. However, we can approximate Levy distributions by introducing a 
truncation at some large value of the argument either in a form of the sharp cut-off or as a rapidly falling tail. Then 
a sum of stochastic variables converges to the normal distribution, according to the central limit theorem. Since 
in the physical phenomena process values are usually finite, introducing truncated distributions is realistic. In the 
random walk theory, the truncated Levy flights are often considered p| |g, [T^-[ l6l . They agree with the Levy flights 
for an arbitrarily large jump value; deviations appear only at very far tails [14j . However, there are also remarkable 
differences between processes which involve the stable distribution and the truncated one. We will demonstrate that 
difference for a simple case of the linear noise. 

Let F(x) = and G(x) — x. If the standard rules of the calculus work - we can expect that for the Stratonovich 
interpretation - the variable in Eq.© can be changed. As a result we obtain from Eq.© the probability density 
distribution in the 'log-Levy' form 

p{x,t) = T-\L a (]n(x/x ),i), (7) 
\x\ 

where L a denotes the Levy distribution with order parameter a, width parameter t and xo = a;(0). If the process 
is continuous the point xq acts as an absorbing barrier, i.e. x > (x < 0) for xq > (xq < 0). It is the case for 
the Wiener process but it may no longer be true if the variance is divergent; then Eq.© is no longer valid. The 
distribution in the form Eq.© for a — 2 is known as the log-normal distribution and it is frequently encountered in 
nature, e.g. electron velocities in the solar wind [l7| . as well as rainfall amounts [H| obey this statistics. Moreover, 
it can serve as a natural model of the multifragmentation [l9j . 

On the other hand, we can solve Eq.© for the Stratonovich interpretation directly from the definition, by means 
of Eq.©. The discretisation gives us x 2 = x\ + (xi + x<i)r]\T X ' a j2\ therefore x 2 = Xi(l + ai)/(l — ai), where 
ai = ?yiT 1 / Q /2. The final solution of the stochastic equation reads 

n , 

x(t)=x l[t±^, (8) 
ti 1 " a * 

where n — t/r. If a cut-off is introduced, ai I for a small r. We take the logarithm of Eq.©, approximate ln(f + x) 
by x and neglect terms of the order r 2 and higher. That procedure yields 

Hx(t)/x Q ) = T^ a ( m + + Vn ). (9) 
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FIG. 1: (Colour online) Probability distribution, calculated from Eq.©, for the case G(x) = x and F(x) = with a = 1.5 
at t = 1, compared with the distribution calculated according to Eq.© with a truncation at a large value of the noise, such 
that Oj > (upper curves for x > 0: solid green line and dashed red line, respectively). Those distributions are identical. The 
distribution which follows from Eq.((5J, but without any truncation, is marked by the dotted black line. The result of the ltd 
interpretation is also shown (dashed blue line). 



The above expression converges to the Levy distribution with the order parameter a, unless x is large compared to the 
cut-off position, and we obtain Eq.(J7J). Since the variance is finite for the process with the truncated distribution, in 
the limit n — > oo the normal distribution must be reached, according to the central limit theorem, but the convergence 
is extremely slow. 

The case of the stable distribution is distinguished by the presence a considerable number of events for which <Zj is 
not small for any given r. Difference in respect to the case of the truncated distribution, due to the presence of those 
events, becomes visible when we consider the distribution p(x, t) for negative x/xq. Obviously in the limit of small r, 
x(t)/xo > if any cut-off is introduced. For the case without truncation, x(t)/xo may turn to the negative. Let us 
estimate the probability P that this cannot happen, i.e. that all factors in Eq.© are positive. We have P — (pi) n , 

where pi — J^ oc p(ai)dai — f_ T p{rf)drj. Since r — > 0, we may insert the asymptotic form of the Levy distribution, 
p{rf) = l/\rj\ 1+a (\r)\ > 1). Then p t = 1 — Jy T i /a pdrj « 1 — r/a2 a . Finally we have 

P=( 1 -- J ^Y -^xp(--M. (10) 



a2 a J \ a2 a 

Probability that at least one of the terms in Eq.® becomes negative, 1 — P, appears finite and it rises with time to 
unity. One can easily demonstrate that P converges to one with r — > for any t if p(rj) is normally distributed. In 
this case p(x, t) = for x < 0. 

Numerical analysis of the above case is presented in Fig.l. The distribution ([7]) was evaluated by means of the 
series expansion 

1 ^r[l + (2»+l)/a] 



L «^ = ri^2^ (2n + l)!! Kt^J ■ 

n— 



The result for the truncated distribution is identical with Eq.([7]) whereas the case without any cut-off (marked by 
dots) exhibits a branch for the negative x. However, both distributions for x > are very similar and then Eq.([7j 
can serve as an approximation of the Levy stable case. The result for the Ito interpretation, Eq.©, is also presented 
in Fig.l. It falls much faster than the Stratonovich one. 



III. NONLINEAR CASE 



In this section we consider stochastic processes which are governed by Eq.((T]) with a nonlinear deterministic force. 
This problem is an important generalisation, compared to the linear case, since the corresponding Newton equation 
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may become nonintegrable and the dynamics is then chaotic. It happens for a periodic time- dependent driving (the 
Duffing oscillator) or if the system has more than two degrees of freedom [2(|. We assume the algebraic F(x) and 
G(x): 

F(x) = -\xpsgnx and G(x) = \x\- e/a . (12) 

In the new variable, 

y(x) = -^U^+^sgnz, (13) 
a + 9 

the Langevin equation, Eq.([S]), has the additive noise. The corresponding fractional Fokker-Planck equation is of the 
form 

^p(V,t) = K^-\yf S gn(y)p(y,t) + ^-p(y,t), (14) 

where fi = 1 — (1 — 7)(1 + 9/a) and K = (1 + 9/a) 13 . The drift term in Eq. (fl"4"|) corresponds to the effective potential 
V{y) ~ |?/|( 7_1 )( 1+8 / Q ). We are interested in the asymptotic shape of a steady-state solution ps(x). The solution for 
large \y\ can be found by taking into account small wave numbers in the Fourier expansion. The Fourier transform of 
Eq. (fT4| in the stationary limit reads 

Kk^{\yf- l p s {y)) = \k\ a p s (k). (15) 



(16) 



We assume the solution in the form of the Fox function [2JJ, [22J , 

(a u ill), (1/2, 1/2) 



Ps (y) = NH& 

(0,1), (6a, 1/2) 



where N is the normalisation constant and the coefficients are to be determined. Some useful properties of the Fox 
functions are presented in Appendix. Eq. (|16j) represents the stable and symmetric Levy distribution for a± = 1 — 1 /a, 
A\ = 1/a and 62 = 1/2 [23J. We insert Eq. pj?)) into Ea.(Ti3|) and apply the general formula (|A1[) in order to get rid of 
the algebraic factor. Then we calculate the Fourier transform, according to the formula (|A2[) . and expand both sides 
of Eq. (fT5|) by using the Fox function series representation, Eq. (|A3j) . Eq. (fl"5|) takes the form 



Kk ~Hk ( ci + c2 ^ wi + c3|fcr + o(fc2) ) = |fc|Q[i + ° (|fc|a)] ' (i7) 

where Wi = (1 — ai)/A\ — fj, W2 = (2 — ai)/A\ — (3 and a are constants. The above equation is satisfied if w\ — and 
W2 = a which conditions determine the coefficients: ai = 1 — /3/a and A\ = 1/a. The condition Kac^ — 1, where 

_ a(a + 1) r(-a)T(a + 8 + l)cos(7rq/2) 
C3_ 2tt r(l + (a-fl)/2)r(-62 + (l-a + fl)/2)' 1 j 

can be satisfied by an appropriate choice of 62. 

The asymptotic behaviour of p s (y) follows from expansion of the Fox function in powers of it can be obtained 

by a variable transformation y y^ 1 by means of Eq. (|A4l) and by applying Eq. (|A3[) . The first term produces 
the result p s (y) ~ |y| _Q_ ' 3 (\y\ — > 00) which, after transformation to the original variable according to the formula 
Ps{x) = p s {y(x))\dy / 'dx\, yields the final result 

p s (x) - \x\-^+ e +^ (\x\ -> 00). (19) 

To satisfy the normalisation condition, we assume a + 9 + 7 > 1 . Eq. (fT9"f predicts the Levy stable distribution with 
a divergent variance for a + 9 + j < 3. 11 a + 9 + j > 3, the variance is finite though higher moments may be 
divergent. Therefore, long tails of the distribution can be confined either by choosing a sufficiently steep potential or 
an appropriate noise. The latter must be such that amplitude declines with position sufficiently fast (large 9) and/or 
the order parameter a is large (steep tails). The case 7 = 1 corresponds to the harmonic oscillator; it is discussed in 
Ref . [l2| . If 9 = and 7 is an odd integer, an analytical expression for the stationary probability distribution, valid 
for arbitrary |x|, can be derived [Hj]. 
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FIG. 2: (Colour online) Stationary probability distributions, for the system given by Eq. (|12[) . calculated by applying the 
transformation (|13[) (lines) and by using Eq. (|20[) (points) for a = 1.8, 7 = 2.5 and two values of 6: -0.2 (the case with a 
maximum in the origin) and 0.5. 



On the other hand, we solve Eq.([T]) by a numerical simulation of stochastic trajectories. It can be performed in two 
ways. First, we directly apply the discretisation formula which follows from Eq.([3]) and is of the form 25] 

x i+ i = Xi + F{ Xi )T + [G( Xi ) + G(.t 1+1 )]t/ 4 t 1 / q /2. (20) 

To find the process value, one has to solve, at each step, the following nonlinear equation 

x i+ i - aiG(x i+ \) - Xi~ F(xi)r - G(xi)a,i = 0, (21) 

where = ^r 1 / Q /2. For that purpose we apply the parabolic interpolation scheme (the Muller method) 26]. The 
algorithm must be carefully implemented since, due to the explicit multiplication of the noise by the x— dependent 
factor, the round-off errors are large and then it is difficult to achieve a high accuracy of the results. Alternatively, 
we can first transform Eq.([T]) to Eq.©, simulate trajectories y(t) to find p(y) and finally transform the distribution: 
p(y) —¥ p(x). Comparison of both algorithms is presented in Fig. 2 for positive and negative 9; the distributions are 
actually identical. In the following simulations we apply the method of variable transformation. 

Distributions which are initially positioned at x — evolve with time to the steady state. The convergence appears 
very fast. Example of the evolution is presented in Fig. 3 for a = 1.5 and 9 = 1. The stationary distribution is reached 
already at t = 1. 

Various sets of parameters a, 9 and 7 define processes which are either stable Levy ones, with divergent variance, or 
processes with heavy tails, for which the variance exists but higher moments are divergent. Examples are presented 
in Fig. 4. The case of negative 9 and a weakly changing F(x) (the case 1. in the figure) corresponds to the slope 2.6, 
in the other cases the slope is larger than three. Slopes of the straight lines in the figure follow from the asymptotic 
formula, Eq. (ll9[) . and they agree with the numerical results. 



IV. ESCAPE FROM A POTENTIAL WELL 



A particular case of Eq.([T]), which involves the nonlinear deterministic force and the boundary conditions, is the 
problem of passing over a potential barrier. This problem is of great physical importance and it has been extensively 
studied for the case of the normal distribution [271 ] . For example, fusion of heavy ions in nuclear physics consists in a 
transfer of mass over the Coulomb barrier. A multiplicative noise emerges when one considers a parametric activation 
of the potential, i.e. if height of the barrier randomly varies [28| . Increasing intensity of the multiplicative noise in 
the bistable stochastic system can produce a stochastic resonance [2!| . Properties of systems driven by general Levy 
stable noises may be different than those for a = 2. In particular, a waiting time for noise-induced jumping between 
metastable states may depend, due to the presence of single long jumps, more on the width than on the height of 
the barrier (30j . Moreover, the ratio of first mean passage time from one to the other minimum is no longer twice 
of the time to reach the top of the barrier [3lj . Asymmetry in the Levy distribution affects the escape time; it can 
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FIG. 3: (Colour online) Time evolution of the probability distribution in the case Eq. (|12[l for a — 1.5, 9 = 1 and 7 = 2.5. The 
distribution was evaluated at the following times: 0.001, 0.01, 0.1, 0.5, and 1, which cases correspond to the rising width. The 
steady state is reached at t = 1. 




FIG. 4: (Colour online) Stationary distributions, calculated numerically, for the system given by Eq. (|12[) (points). The following 
cases are presented: 1. 7 = 1.6, a = 1.5, and 9 = —0.5; 2. 7 = 2.5, a = 0.5, and 9 = 1; 3. 7 = 2.5, a = 1.5, and 9 = —0.5; 4. 
7 = 2.5, a = 1.5, and 9 = 1; 5. 7 = 3.5, a = 1.5, and 9 — 1; (from top to bottom at the right hand side). Slopes of the straight 
lines follow from Eq. lfT^)) . 

both enhance and suppress the escape events. The rate of escape, as a function of the parameter a, is discontinuous. 
It was recently demonstrated that a double stochastic resonance can be observed in a single well potential without 
explicit external driving, if the Levy stable noise is introduced [32j . 

In this section we consider the problem of escaping from the potential well for the multiplicative noise and, in 
particular, the dependence of the mean first passage time (MFPT) on the specific interpretation of the stochastic 
integral, either Stratonovich or Ito. We assume the potential in the form 

V(x) = ^ - fx 2 (22) 

and the noise amplitude G(x) is given by Eq. (jl"2"j) . The main quantity of interest is the dependence of MFPT on the 
parameter 6, which quantifies the noise amplitude variability. The transformed Fokker-Planck equation with additive 
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FIG. 5: (Colour online) The effective potential in the Stratonovich interpretation, calculated from Eq. (|24p . for A = 1, B = 0.1 
and the following values of 9/a: 4/3, 2/3, 0, -1/3 (from top to bottom). 
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FIG. 6: (Colour online) Mean first passage time as a function of 9 for a = 1.5 and 1.2, calculated for the potential 
A = 1 and B — 0.1. Results for both interpretations of the stochastic integral are presented. 



with 



noise is the following 



where 



F(y) = (l + 9/a)y B — A[(l + 6/a)\y\] 



2a/(a+0) 



(23) 



(24) 



The effective drift F(jj) depends only on the ratio 9/a. We can infer some qualitative conclusions about the dynamics 
from the shape of an effective potential, which follows from Eq . (f24| . This potential, as a function of the original 
variable x, is presented in Fig. 5. The height of the barrier falls sharply with 9/a - the potential is very shallow for 
negative 9 - whereas position of the barrier is constant. Therefore we can expect a suppression of transport, which is 
defined by the boundary conditions in the variable x, for large 9/a. 

A numerical analysis of the potential barrier problem must take into account that the system under consideration 
is limited in space. Long tails of the distribution may not manifest themselves if the available space is too small. 
Therefore, in the following, we rescale the system by putting in Eq.@ er = 0.1. MFPT, defined as a time the particle 
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FIG. 7: (Colour online) Mean first passage time in both interpretations of the stochastic integral, as a function of a, for some 
values of 8. 



needs to pass for the first time from the left minimum of the potential to any x > 0, was calculated by numerical 
solving of Eq.© with the initial condition x(0) = —y/B/A and the boundary condition at the absorbing barrier, 
x — 0. The average time, as a function of 8 for two values of a, 1.5 and 1.2, is presented in Fig. 6. It falls sharply with 
8, like an exponential, which is a consequence of the shallow effective potential for large 8. Results for both values of 
a are similar, they differ only by a constant factor. MFPT rises with a since jumps become shorter. This result is 
presented in Fig. 7; a difference between the case of positive 8 and negative one is substantial. 

Predictions for the Ito interpretation are also presented in Figs. 6 and 7. The dependence of MFPT on 8 is much 
weaker than for the Stratonovich case; MFPT falls algebraically for large 9. In this case, the dynamics is affected 
by the noise only near the top of the barrier - the noise is then localised and very strong - while inside the well the 
deterministic trapping dominates. Results for different a are qualitatively the same. Moreover, MFPT rises with a 
in the Ito case, similarly as in the Stratonovich interpretation. 



V. SUMMARY AND CONCLUSIONS 



Levy distribution is characterised by long tails which cause the divergent moments. If noise in the Langevin equation 
is defined in terms of the Levy distribution, the influence of the tails can be confined cither by the potential or by the 
variable noise amplitude. Then the process which is described by the Langevin equation may have a finite variance. 
We studied such processes by solving the Langevin equation with the multiplicative Levy noise and nonlinear drifts. 
The asymptotic shape of the stationary probability distribution depends both on the noise amplitude, assumed in the 
power- law form with the parameter 8, and on the potential slope 7, according to Eq. (|19[) . If the above parameters 
are large enough, the variance is finite for any order parameter a. The asymptotic formula (|19[) is valid only for 
the Stratonovich interpretation; in the Ito interpretation, the distribution is less sensitive on the slope of the noise 
amplitude (the parameter 8). For the case without drift and with the linear drift, the asymptotic formula is the same 
as for the additive noise, i.e. the dependence on 8 does not appear fl2l [l3j. Then the variance is always infinite. 

The difference between both interpretations of the stochastic integral is also visible in the problem of escape from 
the potential well. This problem was studied numerically: MFPT was calculated, as a function of a and 8. The 
effective potential, which includes the variable diffusion coefficient in the Stratonovich interpretation, possesses a 
high barrier when 8 is negative. As a consequence, MFPT rapidly falls with 8. This effect is not observed in the 
Ito case: MFPT falls with 8 according to a power-law. Moreover, MFPT rises with a in both interpretations. The 
above conclusions are valid only if the relative size of the system is large enough to allow the long tails of the Levy 
distribution to manifest themselves. It was ensured in the calculations by taking a small value of the noise parameter 
a. 

Both analytical and numerical calculations for the case of the Stratonovich interpretation were performed by using 
the statement that rules of the ordinary calculus apply and change of variables is possible. That statement is exact if 
the variance is finite, in particular for the truncated distribution. Otherwise, the Langevin equation in the transformed 
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variables, Eq.©, offers only an approximation to Eq.([T]) since for the case of Levy stable processes those equations 
are not strictly equivalent. The approximation is quite accurate but one can also encounter qualitative differences. 
We demonstrated, by considering the case of the linear noise, that the stochastic variable may change its sign, which 
is forbidden for the case of the normal distribution or if the cut-off is present. A possibility to use Eq.© is of great 
practical importance. It enables us not only to perform analytical calculations but also offers a simple numerical tool 
of much higher precision than the direct integration of Eq. ([I} . 



APPENDIX 

In the Appendix, we present properties of the Fox functions which are used in Sec. III. The multiplication rule 





{a p , Ap) 






(a p + aA p ,Ap) 


X 




- H pq 


X 






(b g ,B q ) _ 




(bq+<TB q ,B q ) 



(Al) 



where x > 0, allows us to evaluate products involving algebraic terms. The cosine Fourier transform is given by the 
following expression 



rrmn 





[a p , Ap) 


X 






(b g ,B q ) _ 



cos{kx)dx = -H q ^ p+2 



(l-b q ,B q ), (1,1/2) 
(1,1), (1 - a p , Ap), (1,1/2) 
Numerical values of the Fox function can be obtained by means of the following series expansion 

I\T=i,^ h rfe - UU r(i - aj + Vg) { _ irx ( bh+ ,)/B h 



(A2) 



TTtnn 

pq 





(a p , Ap) 


X 






(b g ,B q ) _ 



EE 



h=1 ,=o IlUra+i r(l - bj + B^) YlU+i r(aj - Ajttf) vlB h 



(A3) 



where n ^ 0. The asymptotic expansion results from the property: 

(a p , A p ) 

Tjnm ' 
= H qp ~ 

{b q ,B q ) ' (l-d p ,Ap) 



I J III II 

H p<i 



(l-b q ,B q ) 



(A4) 
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